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Abstract 



^ ' A general real-time formalism is developed to resum the self-energy 

operator of broken symmetry scalar field theories in form of self-consistent 
t**"* ■ gap equations for the spectral function. The solution of the equations is 

C**) ' approximated with finite lifetime quasi-particles. In the Landau damping 

rates viscosity terms, analogous to gauge theories, appear, what leads to 
a finite damping rate for the long wavelength Goldstone modes. 

ON 

^ ■ 1 Introduction 



Scalar theories appear in many physical situation as part of fundamental the- 
53^1 ories (eg. Higgs sector, inflaton) or as effective models (eg. a - tt system). 

' Understanding their dynamics is a precondition to the clarification of various 

non-equilibrium phenomena in cosmology and heavy ion physics. 

An interesting and probably detectable non-equilibrium process is the forma- 
r> | tion and decay of disoriented chiral condensates (DCC) in heavy ion physics |jj . 

■ The equations of motion which describe the system are studied by many authors 

[^) at classical level and using quantum corrections. The induced current coming 
from the polarization of the thermal medium usually is taken into account at 
one loop level. Higher loops mean incorporation of quasi-particle scattering and 
may have important consequences in certain physical situations. In gauge theo- 
ries, for example, a hierarchically organized two-step one loop analysis resulted 
in the appearance of viscosity terms j|] in the dynamical characterization of the 
lowest energy scales. 

What consequences may we expect from higher loop calculations? Let us 
consider the two diagrams relevant for Goldstone damping shown in Fig. [l]. The 
imaginary part of diagram a) corresponds to the physical process G(p) G(k) — > 
H, where all the particles are on mass shell and G(k) comes from the heat bath. 

*Zoltan Magyary fellow, e-mail: jakovac@planck.phy.bme.hu 
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Figure 1: Cut diagrams contributing to Goldstone damping. G means Gold- 
stone, H means Higgs modes. 



In this process, if p — > then k must go to infinity in order to satisfy the mass- 
shell relations. Diagram b) corresponds to the scattering G(p) H — > G H, where 
H comes from the heat bath. Here, even for p — » all the external momenta 
can be kept finite. In general we will have larger phase space for the processes 
contributing to the damping which will increase the damping rates. This effect 
can be the most dramatic in those cases, where the one loop damping rate is 
small, as in the case of long wavelength Goldstone modes j|, |[ || . 

The other effect we can expect is due to the decay of the intermediate par- 
ticles. At one loop level the intermediate particles are stable, but higher loops 
account also for their own damping, what is another way of interpreting in dia- 
gram b) the inner blob. That means that there is a finite range in time, where 
the interactions are effective. In general this may lead to "loss-of-memory" 
effects for time scales beyond the typical inverse damping rates (cf. also ||). 

In this paper I would like to take into account the dynamical effects of a 
subset of higher loop diagrams by summing up self-consistently the real time 
self-energy functions (Section 2). The resulting generalized gap equations are 
solved in the quasi-particle (Breit-Wigner) approximation (Section 3). The 
results are applied to the calculation of Goldstone and Higgs Landau damping 
and, in particular, the Goldstone on-shell damping rate, at frequencies much 
smaller than the width parameters of the resulting Breit-Wigner shapes (Section 
4) . The conclusions of the investigation are summarized in Section 5 where also 
some important questions left for future studies are discussed. Some formal tool 
used for obtaining results of the main text are presented in an Appendix. 

2 Static and real-time resummations in the spon- 
taneously broken phase of the O(N) model 

The Lagrangian of the model of our investigation is written as 



£($) = -(a$) 2 - - 



(1) 
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where $ represents a vector in the N-dimensional space $ = (4>i,...,$jv)- 
In the symmetry broken phase, where m 2 < 0, the field has a non-vanishing 
expectation value (<f>) = $ =const, which can be chosen to point in the $1 
direction. We introduce a new field </? = $ — $, and the Lagrangian takes the 
form 

C = £($) - C tadpole + ±<pK°(id, m)<p - 9 ^w 2 - §^ 2 ) 2 , (2) 
where the tadpole terms are 

q 2 - 

Ctadpole = m 2 $(fi + ^-$Vi> ( 3 ) 
and the kinetic term is determined by the kernel 

K Q ab (td, ro) = f(id) 2 -m 2 - |U 2 ) 6 ab - ^ 2 6 al 6 bl . (4) 

We will omit from the Lagrangian the purely background dependent and tadpole 
terms in the sequel. 

It is well known that the negative mass squared leads to IR instabilities, 
which appear also in static quantities, for example as a complex free energy 
density for small $'s. In order to avoid these effects one has to resum the per- 
turbation series, which allows for the use of the physical mass in the propagators 
. A consistent method of static resummation consists of introducing thermal 
counterterms 

C = \ V K\id, m T ) V - ^(IvV 2 - f^ 2 ) 2 + \im% - m 2 ) V 2 , (5) 

where the last term is also treated as interaction. The determination of the 
thermal mass can be done with help of gap equations: the masses (for example 
G^ 1 (k — 0)) are stable against loop corrections ||. Since the thermal countert- 
erm preserves O(N) symmetry the Goldstone theorem ensures the masslessness 
of the Goldstone bosons in the broken symmetry case (also at finite temperature 
[|| On the other hand the Higgs mass calculated from the tree level and 
one loop formulae should coincide. These requirements lead to 

a 2 - 

= m 2 + y$ 2 + U l (k =0,m = m T ), 

m 2 " = m 2 + II 1 (fc = 0, to = mj>). (6) 

These are two equations for the two quantities $ and to#. In more complicated 
cases (especially in gauge theories) we can rearrange the perturbation theory in 
a momentum dependent way p"o| . 

In dynamical systems there are more scales than in the static case. In the 
O(N) model besides the temperature (T) and the Debye mass (gT) there appear 
the on-shell damping rates for the Higgs and Goldstone bosons (g 2 T). Especially 



3 



in the Goldstone case, where the Debye scale is missing (ma = 0) can the 
damping play an important role. Also at low external momenta (both ka and 
k falling for instance into the range of Landau damping) these scales may have 
an impact on the results. Therefore it is worth to work out a consistent method 
for the incorporation of the damping scales into the resummation. 

On-shell damping arises when the mass shell falls on some cut of the prop- 
agator. The cuts give momentum dependent contribution to the self-energy, 
which suggest to work with a momentum dependent resummation scheme 

C = ± V K(id)<p - ^(#vV - §^ 2 ) 2 + ±<pP(id)<p, (7) 

where K(id) = K°(id) — P(id). The counterterm has to preserve O(N) sym- 
metry at least for the zero momentum modes in order to be able to satisfy 
the requirements (ip) = and m G = simultaneously to all orders of the re- 
summed perturbation theory. It also has to maintain 0(N — 1) symmetry of 
the Goldstone modes for all momenta. 

At finite temperature in the real time formalism one works with several 
propagators (G < , G > , G c , G a ), and we have to relate them to the kinetic term 
in the Lagrangian. The most comfortable way to do this is first to express 
the spectral function pij — ([<fi(x),<pj(0)]) in momentum space. It is diagonal 
Pijip) = $ijPi(p)i and 

( ) = -2Imifr(p) 
MW (ReK^y + (JmK^y K °> 

where JmK(p) is the shorthand notation for lim £ ^o Im.fr (p + is), and similarly 
for KeK(p). lmK(p) has nonzero values along the cut of the propagator. The 
O(N-l) symmetry of the Goldstone bosons ensures the i-independence of p. 
We can express the propagators through the spectral function as follows 



■ Pun- 
US 



iGf(p) = n(p ) Pl (p) iG>(p) = (1 + n(p )) Pi (p) 

iG^t, p) = Q(t) Pi (t, p) + iGf(t, p) iGf(t, p) = iG> (t, p) - e(t) Pi (t, p). 

(9) 

Now we are in the position to develop perturbation theory with "resummed" 
propagators. The Feynman rules are unchanged (cf. fl]}| ), only the propagators 
are modified. We are interested in the resummed, retarded self-energies (i = 
2, . . . , N and a, b = 1, ... , TV) for which the one-loop contributions can be written 
with help of two types of Feynman integrals (cf. Appendix) 

nf (k) = ^Si + £(N - l)St + ^ 2 Sn(k) + ^ 2 (N - l)Su(k) 

2 2 4 

n f ( k ) = J Si + 9 - (N + l)Si + y* 2 ^(fc). (10) 
The physical masses are 

2 2 

m% = m 2 + ^-$ 2 + n!(0), m 2 G = m 2 + +n t (0). (11) 

2 o 
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Figure 2: Self-energy contributions 



The quantities S a and S ab correspond to the diagrams shown in Fig. ^|. S a 
is momentum independent, it yields merely a modification to the mass. This 
contribution is proportional to G c (x = 0) = G K {x = 0), since p(x = 0) = 0. In 
Fourier space it reads 

d 4 p 

The other diagram is nonlocal, especially it has imaginary part in the momentum 
representation. In the coordinate space 



S a := I 7 ^n( P0 )p a (p). (12) 



S ab (x) = (iG c a (x))(iGt(x)) - (iG<(x))(iG<(x)) 

= Q{t)[p a (x)p b (x) + Pa {x)iG<{x) + Pb (x)iG<(x)}. (13) 

In Fourier space it has the form 

o ru\ f d Po d Po d3 P Pa(PQ'P)Pb{Po,k-p) , 

Sab{k)= , - 3 — ——(l + n a (po)+n b (p )). 14 

J 2n 2-k (2-kY k ~p -p' + ie 

Its imaginary part with respect to the kg variable is given by 

ImS ab (k) = -\ / ( 2^4 (2^)4 Pa( P )p b (p')(l + n a ( Po ) + n b (p' ))(27r) i 6(k-p- P '). 

(15) 

Using 

1 + n a (p a ) + n b (p Q ) = (1 + n a (p Q ))(l + n b (p' Q )) - n a (p Q )n b (p' ) 

the above result can be interpreted as the difference of pair creation rates into 
(and annihilation out of) the states a and b where the density of states are p a 
and p b @. 

The expectation value for the field operator can be worked out very simply. 
(<Pi) = is trivially fulfilled, (ipi) = requires 



$[ m -» + *-& + ->-S 1 + ^(N-l)S i ) =0. (16) 



Using this relation with m G — leads to 



G 

S 1 + S l + ^ 2 S ll {k = 0) = 0. (17) 
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This is fulfilled in the perturbation theory H, but in a general resummation 
might fail. In order to maintain the Ward identities we have to make some 
further efforts, and perform also vertex resummation |l3|. In this context it is 
sufficient to add and subtract from the original Lagrangian a piece of the form 
A<? 2 <I>f/24, where Ag 2 ~ g 4 . The original Lagrangian does not change, but, if 
we take into account the two pieces at different loop orders of the perturbation 
theory, their contributions break the O(N) symmetry at each level. At tree level 
the Goldstone mass is not affected, but in the second term of (|l6| ) leads to the 
substitution g 2 <& 2 — > (g 2 + Ag 2 )<& 2 . At one loop the compatibility of (tp) = 
and m,Q — is fulfilled if 

An 2 - a 2 - 

-y&-S 1 + Si + Z-* 2 Su{k = 0) = 0, (18) 
2g z 3 

which determines Ag 2 . Keeping this in mind we omit vertex resummation, since 
it has no effect up to the order of the perturbation theory which is considered. 
Introducing the notations 

2 2 

m 2 T = m 2 + L.S 1 + L.(N-l)S i , 
2 6 

9t = 9 2 (l + yS n (fc - 0) + ^(N - l)S u (k = 0)) (19) 

(<7y, as can be easily checked, is the one-loop corrected coupling) we have the 
relations 



9 3 g 2 



Therefore to the given order we can substitute 5 4 $ 2 by Sg 2 m 2 H . 

For the optimization of the time-dependent perturbation theory we can use 
the same ideas as in the static thermal counterterm case: we demand that 
the one loop propagator be identical with the tree level one for all momenta. 
Implementing this requirement in the Schwinger-Dyson form for the inverse 
propagator 

G-\ P ) = K{p) - n(p;K) + P(p) = K{p) (21) 

yields 

Pip) = Il(p;K), (22) 

ie. the dynamical counterterm is the self-energy calculated using the one-loop 
exact propagator. Earlier we have determined the expressions for the retarded 
self-energies as functionals of the spectral functions, and also the spectral func- 
tion is expressed through the kinetic term K. Thus we can write down a com- 
plete set of renormalized functional equations 

Hi(fc) = ^Si + £(N - l)Si + ^-m 2 H S 11 {k) + ^m 2 H {N - l)S u (k) 

2 2 2 

Hi(k) = ^-S 1 + ?-(N + 1)5, + L.m 2 H Su(k) 
6 6 3 
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-2Imn 1 (fc) 

Pi(k) 



(k 2 -m 2 H - Re(n x (fc) - IIi(O))) 2 + (Imlli(fc)) 2 
(k] = -2InJL;(fc) 

where mu comes from eq. (|20|). This set of functional equations can be, in prin- 
ciple, solved. Their construction represents the main conceptual result of our 
investigation. Now, we shall proceed with developing practical approximation 
schemes applicable to these equations. 



3 The Breit-Wigner approximation 



We can simplify the solution of our gap equations in (|23|) with some assump- 
tions. The most severe is the quasi-particle approximation, that is we assume 
that numerically the most important contribution comes from a narrow region 
around the mass shell w(k) of the spectral function (cf. Jljl). This assump- 
tion means that we take into account only those effects which stem from the 
resonant broadening of the mass-shell. If contributions coming from the mo- 
mentum ranges of the spectral function which are far from the single particle 
mass shell prove to be important, then our Breit-Wigner approximation is not 
valid any more. A useful characteristics of the actual situation is the value of 
the wave function renormalization constant Jig] . If its value is far from 1, 
the quasi-particle approximation fails, while for Z « 1 we can hope that we are 
on the right track. Formally this assumption corresponds to the Breit-Wigner 
approximation to Ii{k), with ImiT(A:) « — 2wk7k independent on kg. For 7 <C lo 
we can write 

Pa W (p) = 7~T ( 5 7a( P )(Fo -w (p)) - <5 7a(p )Oo +u> a (p))) , (24) 

where w a (p) will be approximated by the mass-shell relation as 

u, 2 (p)«p 2 +m 2 , (25) 



m 2 standing for the resummed masses (cf. @), ie. we neglect the momentum 
dependent real part of the self-energies. Furthermore, 

8 7 (x) = -lm^— = --^-^ (26) 

is a smeared out delta function, where the on-shell damping rates 7i(p) yield 
the width of the functions. In the free case we get back to the usual expression 

lim pf w (p) = (2n)e(p Q )5(p 2 m 2 ). (27) 

7-»0+ 

Under these approximations we can simplify the gap equations using instead 
of the last two equations of ( ^3|) eq. (0) with the usual form of the on-shell 
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damping 



7i (P) = 



-Imnf(o;i(p),p) 
2wi(p) 



7i(p) 



-Imnf(|p|,p) 
2|P| 



(28) 



As a further approximation we treat 5 7 acting on any function smooth 
around the mass shell as an ordinary Dirac-delta distribution. This is appro- 
priate if the rate of variation of the smooth function is slower than (or at least 
comparable with) 7. Since 7 is the smallest scale in the system (7 ~ g 2 T) 
we can apply it in many cases. If we substitute the Bose distribution with its 
on-shell value, S a turns out to be 7-independent 



P — [l + 2n(w« ]. 



(2tt)3 2uj a 

For the computation of S a b we introduce the notations 

p' = k-p, 7 = 7 a (p), 7' = 7b(p'), w = w a (p), u'=u)btj?') 
n = n(uj), n' = n(uj'), 5^ = 6 a {p ±uj) 5 b = 5 b (p' ± u'), 

which allows us to write (cf. ( |l4| ) and (p4)) 



(29) 



(30) 



dpo rfpo ^ 3 P 71-2 
"27r"27T(27r) 3 wu7 



(l + n + n ) 



+ t «- j -- (n'-n) 



(31) 



Now we can perform the po and p' integrals with help of the formula (cf. Ap- 
pendix) 



1 

27 



dx 



(x + ia)(Q — x + £7) 



■(7)6(70:) 



fi + i(o + 7 )' 



(32) 



This yields 



dp 



fl — Po + ie 



dpo 



1 



1 



27TZ \Po ± 
1 
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2 f2±w + i(e + 7)' 



Po ± w + «7 / £1 - pa + ie 
(33) 



ie. S a behaves almost like an ordinary Dirac-delta in this case, but its width 
is added to the imaginary part of the denominator. The same formula can be 
applied again to the p' a integration, and finally we arrive at 



S a b{k) 



1 



l + n + n' 



l + n + n' 



(2tt) 3 Aujuj' [ k - uj - uj' + iT k Q + uj + uj' + iT 
n — n' n — n' 



ka + uj — uj' + ir fco — uj + uj' + iT 



(34) 
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where r = 7 + 7'. 

The result formally corresponds to changing the imaginary piece is in the 
perturbative denominators to iT, which has the physical content of replacing 
the adiabatic switching off the interactions by the effect of the particle decay. 
This leads to the realization of our conjecture in the Introduction, namely the 
loss-of-memory effects. After a time ~ 1/T the interactions become ineffective 
and we can expect in the time evolution of any physical quantity the appearance 
of effectively localized (instead of long tail) kernels. 



4 Imaginary part of the self-energies 



The other expected effect, the larger available phase space of final states con- 
tributing to quasi-particle damping will be demonstrated via the actual calcu- 
lation of the imaginary parts. From eq. (23) we conclude 



Imnf (fc) 



2 

„2 



m 2 H ImS u {k) + ?-m 2 H {N - l)Im^(fc) 



Then eq. 
lmS a b(k) 



Imnf (fc) = ^m%lmS u {k). 

UJ) gives 

d 3 p 



(2tt) 3 Aujuj 1 



(1 + n + n')(5 T (k - oj - oj') 
+ (n - n')(5 r (k + uj - uj') - 



- 6 r (k 
S r {k - 



(35) 



V)) 



W)) 
(36) 

We can use again our assumption that Sr, acting on a smooth function be- 
haves as an ordinary Dirac-delta function. Here we apply it to Bose-Einstein 
distributions as well as to damping rates. The latter are functions of \p\, or u> a , 
thus we replace the argument of 7(0/) by the value dictated by the Dirac-delta 
condition. We use also the identity 1 + n(x) + n(—x) = 0, and we obtain 



lmSab{k) = - 1 "wf 4^7 1 



(I + n(u>) + n(ko — u)))Sr(ko — lu — uj) 



- (n(ui) - n(k + u))5r{k + uj + uj') 
+ (n(uj) - n(k + oj))5 r {k + uj - uj') 
+ (n(uj - fc ) - n(u)))Sr(k - uj + uj')] ■ (37) 

The integrand depends only on uj and a/, therefore it is worth to change the 
integration measure correspondingly. In spherical polar coordinates d 3 p = 
2irdpp 2 dx, where x — pk. We change variables from p to uj and from x to 
uj' . Then we obtain in a straightforward procedure 

00 u+ 

d 3 p 



[2irY Aujuj 



7/(^0 



I 

167r|k| 



duj I duj' f(uJ, uj'), 



(38) 
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where 

Lu 2 ± = (p±\k\) 2 +m 2 , u 2 =p 2 +m 2 a . (39) 

The sole (//-dependence remained in the <5r functions. Introducing smeared 
theta functions via the relations 

X 

11 X 

dy S r (y) = Or (x) = - + - arctan - (40) 

Z 7T 1 



(lim r ^ @r{%) = ®(a0) an d 

If (x — a) (x — b) 

9r(a < x < b) = 0r(x — a) — 0r(a; — &) = — I arctan — — arctan 



we immediately find 



7T v r r 

(41) 



lmS ab (k) = — — ,. ; [ dw\ (1 + n(u>) + n(k - w))9 r (w- < k - lu < lu + ) 
167r|k| 7 

— (n(w) — n(fco + w))0r(w_ < — fco — w < w+) 
+ (n(w) — n(ko + u))Qr(u>- < kg + u> < lu + ) 
+ (n(ui — fco) — n(w))0r(w- < w — fco < <•*->+)] . 

(42) 

This shows the other expected effect conjectured in the Introduction, the en- 
larged phase space contributing to the damping, since 0r is a smeared version 
of 0. Next we present some direct consequences of the application of eq.(ff2|). 

4.1 On-shell damping rate of Higgs bosons 

If the coefficient of g 2 T in the perturbative expression of the damping rate is 
0(1), then the corrections from the finite width are negligible, or more precisely 
they cannot be assessed reliably from the present approximation. In these cases 
the Higgs mass is the relevant scale which determines the damping rates, the 
finite width plays only subleading role. As an example let us take the on-shell 
damping of the Higgs component. This phenomenon is the result of the decay 
of the Higgs boson into two Goldstone bosons in the zero width case. Now we 
want take into account the effect of broadening of the resonances in this process. 
Using (||), @ and the first line of @ we write (k 2 = k 2 + m 2 H ) 

oo 

9 2m2 H{N — 1) f . 

l\on- shdi k = inQ , I, I / dp (1 + u p + n fco - P)) 
1987rfc k J 



xe r (|p-|k|| <fc -p<p+|k|).(43) 
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If |k| — » (and p > |k|) then using eq. (|l|) 

lim TT-j-eriw- < s < w+) = 2 — = — d r (^b - s). 

|k|->o |k| dp diVb uj b 

Now uib = p and s — m — p, therefore 

oo 

g 2 m H (N- 1 



71, 



Ml(k - 0) = g mg Q f = 1} (1 + 2»&) / dp6 r (m H 2p). 
9o7r 2 7 



(44) 



(45) 



The integral is formally (1 — r(mjy)/m#)/2, however, within our approximation 
scheme it can be also be approximated by 1/2 (by the rule of the application of 6r 
to a smooth function) . The difference can be tracked back to the neighborhood 
of the p = lower limit of the integration, which yields 5r(mH)- After closer 
inspection of the two-particle cut contributions to this process, it can be seen 
that the far off shell contributions (for example the ones coming from Sn) are of 
the same order of magnitude. As a consequence we cannot give the correction to 
the one-loop result for ^fx,on— shell 

(k), as it was stated above. The conclusion is 
that better approximation scheme is needed to the solution of the self-consistent 
equations in order to improve it beyond the one-loop result in case of Higgs 
bosons. 



4.2 Landau damping and on shell damping rate of Gold- 
stone bosons 

The situation is different in case of the on-shell damping of Goldstone particles 
at small momenta. From perturbation theory Q 



1 / m 



pert 



> 2 H 



which is exponentially small. Here it makes sense to examine, how its self-energy 
is modified due to the quasi-particle width. In our approximation, using (p8[), 
~35|) and the Landau damping part of (M2]) 



oo 

Imn <<*> —SsiS/* 



{n{p) - n(p + fco))@r(w_ < p + k Q < lu + ) 



+ (n(p - k ) - n(p))&r(u>- < p - k a < u> + ) 
If |k| <C rrih, T we can use eq. and write (to 2 = k 2 + m 2 H ) 

oo 

g 2 m 2 H f dpp 



(46) 



Imn^yfc) 



24tt J uj 
o 



(n(p) - n(p + fco))£r(w —p — ko) 



+ (n(p - k ) - n(p))6 T (uJ -p + k ) 



(47) 
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If also fco <C ma, T is fulfilled we expand the Bose-Einstein factors in powers of 
their arguments, and find for the leading term 



t tt ru\ 9 2m H k Q fdpp , dn \ x , \ r,c,\ 
lmU t (k) = -f— / (- — )5r(w-p). (48) 

1/7T J U! dp 



For r = the delta function can never be satisfied. For large p, however, 
uj — p pa m 2 H /(2p), the argument of Sp approaches zero (unlike the previous 
cases, where the argument was constant). Moreover, for large p, T ~ — ImII/2p, 
remains smaller than u> — p. Assuming that this region dominates the integral 
(to be confirmed a posteriori) we neglect T in the denominator and write 

oo 

o 

The integrand receives contribution from momenta p ~ T > niu, therefore the 
starting assumption was correct. Approximating T by a constant, the integral 
can be performed 

ImrF(fc) = - g2 J 2 2° r = -Vi ko- (50) 

Here T = 71 (p) + jj (p) must be taken at the dominating value of the momentum 
p T. Our notation emphasizes that a viscosity term has appeared, similar 
to the one observed in gauge theories || . According to our result the viscosity 
coefficient r] is proportional to T(T), the quasi-particle damping about T. The 
above formula also says that long wavelength Goldstone modes are damped in 
contrary to the zero width result, and 

1i,on-shell{k = 0) = y. (51) 

If we take for T only the Higgs contribution, as calculated from perturbation 
theory § 

g^TlN — 1) T 

4327T TJlH 

4.3 Higgs Landau damping 

Landau damping of the Higgs fluctuation modes in the usual perturbation theory 
comes from the k 2 < contributions of Sn and Su. We now examine their 
modification due to the finite lifetime of quasi-particles. 

First let us consider Su. We will apply high temperature expansion on the 
results. More precisely in the last two lines in eq. ( fi"2"|) we assume p»m, |k|, 
thus neglect masses and uj± as p ± |k|. Thus 

e r (w_ < w ± fc < w+) pa e r (fc + |k|) - e r (fc - W) = e r (-|k| <k < |k|), 

(53) 
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depends on lu only through L. If T = this is nonzero only for fco < |k| and 
the usual abrupt end of the Landau damping phenomenon occurs for fco = |k|. 
For finite width, however, the range of effect in fco might broaden. The Landau 
damping part of ImSii then reads 

oc 

lmS Land.dam V . ^ J dio (n(w - k ) - n(cu + fc ))O r (-|k| < fco < |k|). 

(54) 

The smooth difference n(cu — fco) — n{uj + fco) ~ ~2kodn/duj enhances contri- 
butions from small momenta. Therefore if we approximate T(p) = 271 (p) by a 
constant, we have to take it at a typical momentum p ss 0. From perturbation 
theory (see eq. (ff5|)) 

F«2 7l (0)«g. (55) 
Then high temperature expansion yields 

lmS Land. damp . w ln mg + ^° e r (-|k| < fc < |k|) 

167r|k| m H — fco 

' /C ° r e r (-|k| < fc < |k|). (56) 



8nmH |k| 

This is the zero width result simply modified by the substitution 6(|k| — fco) ~ * 

e r (-|k| < fc < |k|). 

In case of the Landau damping contribution from ImSa the situation is 
complicated by IR divergences. To be on the safe side we introduce there an IR 
cutoff |k| < A < n%H, and use cut Bose-Einstein distributions || 

h(uj) = Q(lu — A)n(w). 

The result of the calculation is thereafter interpreted as induced self-energies for 
an effective model for modes below the cutoff ||, |[ . Then the Goldstone modes 
do not contribute, and the Landau damping part of the Higgs self-energy reads 
(cf. @) 

ImU R Mdamp{k) = „^_ T m H ^6 r (-|k| < fc < |k|). (57) 
167r |k| 

Qualitatively new phenomenon appears in comparison to the zero width case if 
fco, |k| < T, ie. below the scale ~ g 2 T. There (cf. ©) 9 r (-|k| < fc < |k|) w 
2|k|/(7rr), and 



Im nf.^^(| k |, fc < r) = -|4 ^ k = -mfco. (58) 

8tt z 1 



This is independent of k (also the usual kinematical restriction for the Landau 
damping, fco < |k| does not hold), and is valid till fco < T. If we consider it 
as part of the effective action for the low frequency modes, in real space this 
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turns into an ~ idt viscosity term. The presence of a similar effect has been 
demonstrated in gauge theories ||. The present calculation shows that this 
phenomenon is common in effective dynamical theories of the low frequency 
modes of any non-linear field theory. 

5 Conclusions 

In this paper we made an attempt to incorporate higher loop effects into the 
perturbation theory and work out some of the most dramatic consequences they 
might provide. We expect, on general grounds, two types of effects. The first is 
the result of the larger available phase space for the processes contributing to 
the damping which increase the physical damping rates. The other effect is due 
to the instability of intermediate particles which lead to loss-of memory effects. 

We have demonstrated, how to resum self-consistently a subset of Feynman 
diagrams (self-energy corrections), what leads to the closed set of functional 
equations (gap equations) of eq. (|2^). Its general solution is very complicated, 
therefore we have made some simplifying assumptions. Among them the most 
important was the Breit-Wigner quasi-particle approximation. Then the for- 
malism is similar to the zero-width calculations, just the Dirac-delta and the 
Theta distributions are replaced by smeared versions. We have calculated the 
imaginary part of the self-energies in this approximation. 

We have shown, how the Landau damping regimes of the different excitations 
are modified due to the finite width effects. In general we have arrived at a 
formula for fc , |k| <C 7 

Imn Landau damping (fc) = (5Q) 

which corresponds to r\ a dt in the real space. Such a viscosity term appeared 
also in the gauge theories ||. This again confirms our expectation about the 
loss-of-memory effects beyond the time scale ~ 1 /j. Long wavelength Goldstone 
bosons (fc = |k| <C 7) fall into this regime, thus having a finite damping rate 
j° n ~ sheU (k) = T)i/2, where % ~ ji(T), the on-shell damping rate of the Higgs 
bosons with momentum about T. Other on-shcll damping rates which in the 
perturbation theory are proportional to g 2 T with a coefficient of 0(1) receive 
0(7/771) ~ g correction, which is, however, beyond the scope of the present 
approximation. 

After this calculation one faces the question still awaiting clarification, which 
subset of Feynman diagrams yields the correct time evolution of the physical 
quantities. We know, for example, that the kg = Goldstone mode which 
corresponds to the homogeneous rotation of the vacuum is not damped (cf. 

^)) m contradiction with the present result. Probably this puzzle can be 
resolved by a resummation of a subset of vertex diagrams, similarly to the 
resummation in the $ 4 theory ]T6| which leads to an important modification in 
the viscosity coefficients. It is also worth to study, how the present ideas can be 
applied to gauge theories and to systems far from equilibrium. 
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Appendix 



The appearance of retarded self-energy can be determined from the Schwinger- 
Dyson equations for the matrix propagator. We will use here the matrix nota- 
tions, G 11 = G c , G 22 = G\ G 12 = G< and G 21 = G> 



G ij = G o,ij + G o^k u ki G h, 



(60) 



where the indices can be 1 or 2. Writing out the components of these equations 
one finds 



Qil = Q0.il + QO.iljjIj Qjl + G 0,i2 n 2j Qjl 
Qi2 = G Q,t2 + QO.iljjIj Qjl + G 0,i2 n 2j G j2_ 



(61) 



Subtracting the two equations and using G R = G jl — G j2 for both j — 1, 2 we 
find 

G R = G 0.R + fjCil pll + n l 2) G R + G 0, i2(n 21 + n 2 2) G R (g2) 

Subtracting again the equaitons i = 1 and i = 2, then multiplying the result by 
(Gh)- 1 and using G 11 - G 21 = G 12 - G 22 follows 



n 11 + n 



12 



n 21 + n 



22 



0. 



Introducing H R = IT 11 + II 12 we arrive finally at 

G R = G 0,R + G 0,R U R G R 



or 



0,il\-l 



(G ) = (G u - K ) 



(63) 

(64) 
(65) 



Schwinger-Dyson equations. 

The formula under ( |32| ) can be derived in the following way. For 7 > 



1 

2^ 



fix 



(x + ia)(Q — x + £7) 



2tt n- 


- i(a - 


h7) 


1 


1 




2tt fi- 


- i(a - 


r-7) 


1 


1 





fix 



1 



1 



27r £1 + j(a + 7) 



a; + ia. x — il — 27 

oo+ja 00—27 



fix 



fix 



[ln(-l - ia0+) - ln(-l + i 7 0+)] = 



6(a) 



+ j(a + 7)' 



(66) 
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On the other hand 

oo oo 

dx If dx 



2t: J (x + ia) (£1 — x — ij) 2n J (x — id) (— f2 — x + vf) 



®(- a )-^ T \' ( 67 ) 



The two equations lead to the result of eq 
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